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ABSTRACT 

We present results from a detailed spectrophotometric analysis of the blue compact dwarf galaxy 
Mrk 35 (Haro 3), based on deep optical (BVRI) and near-IR (JHK) imaging, Ha narrow-band obser- 
vations and long-slit spectroscopy. The optical emission of the galaxy is dominated by a central young 
starburst, with a bar-like shape, while an underlying component of stars, with elliptical isophotes and 
red colors, extends more than 4 kpc from the galaxy center. High resolution Ha and color maps allow 
us to identify the star-forming regions, to spatially discriminate them from the older stars, and to 
recognize several dust patches. We derive colors and Ho parameters for all the identified star-forming 
knots. Observables derived for each knot are corrected for the contribution of the underlying older 
stellar population, the contribution by emission lines, and from interstellar extinction, and compared 
with evolutionary synthesis models. We find that the contributions of these three factors are by no 
means negligible and that they significantly vary across the galaxy. Therefore, careful quantification 
and subtraction of emission lines, galaxy host contribution, and interstellar reddening at every galaxy 
position, are essential to derive the properties of the young stars in BCDs. We find that we can 
reproduce the colors of all the knots with an instantaneous burst of star formation and the Salpeter 
initial mass function with an upper mass limit of 100 A/ Q . In all cases the knots are just a few Myr 
old. The underlying population of stars has colors consistent with being several Gyr old. 
Subject headings: galaxies: dwarf - galaxies: evolution - galaxies: individual (Mrk 35) - galaxies: 
starburst - galaxies: stellar content 



1. INTRODUCTION 

Blue Compact Dwarf (BCD) galaxies are narrow 
emission-line dwarfs that are undergoing violent bursts 
of star- formation (Sargent & Searle 1970). They are 
compact and low-luminosity objects (starburst diameter 
< 1 kpc; Mb > —18 mag), often with low-metal content 
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(Z e /50 < Z < Zq/2) and high star-forming (SF) rates, 
able to exhaust their gas on a time scale much shorter 
than the age of the Universe. Initially it was hypothe- 
sized that BCDs were tr uly young galaxies, forming their 
first generation of stars (Sargent & Searle|1970 Lequeux 



fc ViallefondTT980l |Kunth k Sargent|l986| |kunth, Mau 
rogordato, fc Vigroux||1988p , but the subsequent detec 
tion of an extended, redder stellar host in the vast ma- 
jority of them has shown that most BCDs are old sys- 
tems (Loose fc Thuan|1986l|Telles|1995||Papaderos etal 
1996a[ |Cair6sP000| | Cairos et al.|fi001a[ =C0 1a; |Cair6s 
et al. 12001 b| = C01b;TCaIr6s et aU2002| =C02; |Cair6s et 



al.|2003[ =C03; |Bergvall fc Ostlin|2002| >. Although so far 
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the SF in BCDs has been commonly described as a burst- 
ing scenario — short int ense episodes of SF separated by 
long i nactivity periods (Thuan 1991 Mas-Hesse & Kunth 
1999]) — , there is increasing evidence that BCDs could 



rather have a gasping star formation, with long episodes 
of activity, separated by short quiescent intervals if any 



2005). 



(Legrand 2000 Legrand et al. 2000 Recchi fc Hensler 



The evolutionary status of BCDs and their star forma- 
tion history, as well as the mechanisms that trigger the 
recurrent star formation episodes, are still open ques- 
tions, and, although much work has been done in the 
last years in this topic, no conclusive results have been 
achieved. 

A fundamental step to approach such questions is a 
detailed analysis of individual nearby objects, using high 
quality data. In fact, most of the work done so far has 
focused on statistical analyses of BCD samples, and only 
recently a few studies have been devoted to examining 
in detail the c h aracteristics of individual objects. |Pa-| 
paderos et al.| (|1998[ |1999| carried out a spectropho- 

335-052 and Tololo 65, two 
BCDs belonging to the iO class (extremely compact and 
low-metallicity objects with no eviden ce for a substan- 
tial u nderlying stellar component; see |Loose fc Thuan 
1986 for the classification scheme); |Noeske et al. " p000| ) 
studied the i,IC galaxies Mrk 59 and Mrk 71 (" Cometary 
BCDs" , galaxies with star formation concentrated on one 
side); Tololo 1214-277 and SBS 0940+544 (als o classified 

T 



as i ,IC) w ere studied by Fricke et al. (|2001 



al 
wit 



2001 ), respectively. Despite the fact t 



H) and lGus 
that nE (o 



Guseva et 
bjects 



a single central starburst superposed on a regular 
redder outer envelope) and iE (objects with a complex 
inner structure consisting of many SF knots atop a regu- 
lar redder outer envelope) BCDs are the most common, 
there are very few studies of similar quality for galaxies 
belonging to these groups, the most notable one being 
the comprehensive an alysis of the properties o f th e iE 
Mrk 86 carried out by [Gil de Paz et al.| (|2000| and [Gil 
de Paz, Zamorano, fc Gallego| ( |2000p 



i'his prompted us to start an extensive observational 
project, whose aim is the thorough analysis of a sample of 
nearby BCD galaxies selected among those already stud- 
ied by our group (COla and COlb). In a first paper (C02) 
we introduced the method used to discriminate and ana- 
lyze the stellar populations in BCDs, and presented the 
results for the iE galaxy Mrk 370. Here we apply the 
same method to Mrk 35, a BCD also belonging to the iE 
morphological class. 

Mrk 35 has an absolute magnitude Mb = —17.75 
(COlb), and is located at a distance of 15.6 Mpc 1 . Ac- 
cording to the morphology and distribution of its SF 
knots, Mrk 35 belongs to the Type hi (Chain/Aligned 
starburst) group, as introduced in COlb. The basic data 
of Mrk 35 are shown in Table [jj . 

Deep optical surface photometry, in the optical and in 
the NIR, was presented in COlab and in C03. The galaxy 
displays the typical surface brightness profile (SBP) of 
BCD galaxies: at high and intermediate intensity lev- 

This distance has been computed assuming a Hubble flow, 
with a Hubble constant Ho = 75 km sec -1 Mpc -1 and applying 
the correction for the Local group infall into Virgo. 

2 A collection of color maps of this galaxy can be found at: 
|http : //www . iac . es/p royect/GEFE/BCDs/ 



TABLE 1 

Collection of data on Mrk 35 



Parameter 


Data 


Reference 


Other names 


Haro 3, NGC 3353 




R.A. (J2000) 


10 45 22 




Decl. (J2000) 


55 57 37 




D (Mpc) 


15.6 


(1) 


A B 


0.031 


(2) 


M B 


-17.75 




m B 


13.21 


(3) 


ray 


12.60 


(3) 


m R 


12.37 


(3) 


mj 


11.84 


(3) 


mj 


11.31 


('D 


m H 


10.83 


('!) 


m K, 


10.43 


('1) 


M HI 


0.46X10 9 


(5) 


M T 


1.5 x 10 9 


(5) 



Note. — (1) Distance computed assuming a 
Hubble flow, with a Hubble constant Ho = 75 km 
sec -1 Mpc -1 , and applying the correction for the 
Local Group infall into Virg o. (2) Absorption co- 
efficient in the B band, from Schlcgcl, Fink b"emer,| 
|fe Davis] | |1998} . (3) Asymptotic photometry ob- 
tained by extrapolating the growt h curves, and 
correct ed for Galactic extinction | |Cair6s et al.| 
|2001b[ l. (4) Isophotal magnitudes (within the 
isophote at 23.0 mag arcsec - 2 in J, 22.0 in H 
and 21.0 in K s ; C03). (5) Neutral hydrogen mass 
Mm and total mass Mt in u nits of Mq; both 
from |Thuan k. Martin| ( |1Q81| ) ■ Units of right 
ascension are hours, minutes, and seconds, and 
units of declination are degrees, arcminutes, and 
arcseconds. 

els the profile is dominated by the starburst component, 
which has quite a complex shape, whereas in the outer 
parts the SBP traces the luminosity structure of the 
old stellar populat ion. Mrk 35 is included in the 



Maz- 



zarella fc Borosoh] ( 1993 1 sample of Markarian ga 
with multiple nuclei. The presence o f Wolf-Rayet stars 
( Steel et al.||1996 Huang et al.||1999 l indicates that the 
galaxy underwent a starburst episode within the last 3-6 
Myr. The caus e of this star - forma tion episode is un- 
clear. Whereas Steel et al. (19961 conclude that the 
star -form ation is likely auto-induced, |Sanchez-Portal et| 
[ah] ( |2000[ ) claim that Mrk 35 show s clear signatur es of a 
merger episode. A recent paper by Johnson et al. (20041 
studies the central starburst by means of near infrared 
and radio observations, and concludes that the hypothe- 
sis of a small scale interaction can not be ruled out. 



2. OBSERVATIONS AND DATA REDUCTION 
2.1. Long-slit spectra 

Long-slit spectra of Mrk 35 were obtained in Febru- 
ary 2002 at the Observatorio del Roque de Los Mucha- 
chos (ORM) on the island of La Palma, using the 4.2m 
William Herschel Telescope (WHT). Observations were 
made with the blue arm of the ISIS double beam spec- 
trograph, which was equipped with a 300 groove mm -1 
grating and a CCD array of 2100 x 4200 13.5/^ pixels, 
a combination which gives a linear dispersion of 0.86 A 
per pixel, and a spectral range of 3600 to 6920 A. A slit 
4 arcmin long and l'/2 wide was used. The positions of 
the slit we used is shown in Fig. [T] over-plotted on a 
continuum-subtracted Ha map. 

The exact positions of the slit were derived from a di- 
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rect image obtained from the ISIS acquisition TV. The 
position angle of the two slit positions was set to 42° 
and 38° (SI and S2 in Fig. [l]and hereafter), so as to pass 
through the brightest central knots and the knots in the 
east tail. We took three exposures of 1800 seconds in the 

51 position, and two exposures of 1800 seconds in the 

52 position. The average seeing during the observations 
was about l'/2. 

Data reduction was performed using IRAF standard 
tasks: after subtracting the bias and flat-fielding the 
spectra, they were calibrated in wavelength. The sky 
spectrum was derived by averaging the signal in two win- 
dows 60-70 pixel wide outside the region where the ob- 
ject is still detectable, and subtracted out. The spectra 
were then corrected for atmospheric extinction and flux 
calibrated by means of observations of spectrophotomet- 
ric standards. 

2.2. Hot imaging 

Images centered on the Ha line and adjacent contin- 
uum were acquired in March 2002 at the Nordic Opti- 
cal Telescope (NOT; ORM), using ALFOSC (Andalucia 
Faint Object Spectrograph and Camera). The detector 
was a 2k x 2k, 15/xm pixel Ford-Loral CCD, which, with 
a scale of 0.188 arcsec per pixel, provides a field of view 
of 6'A x 6'A. The seeing varied around 1'.'2. The complete 
log of the observations is shown in Table [2j 

The image reduction was carried out using IRAF. Each 
image was corrected for bias, using an average bias frame, 
and was flattened by division by a mean twilight Hatfield 
image. The average sky level was estimated by com- 
puting the mean value within several boxes surround- 
ing the object, and subtracted out as a constant. The 
frames were then registered (for each filter we took a 
set of three dithered exposures) and combined to obtain 
the final frame, with cosmic ray events removed and bad 
pixels cleaned out. 

The Ha frame contains the emission from both the Ha 
emission line and the underlying continuum; in order to 
obtain the net Ha we need to subtract the continuum flux 
(the off-line image) from the Ha total emission. First, the 
sky-subtracted off-line frame was rcscalcd in intensity to 
match the continuum in the sky-subtracted Ha frame. 
The scale factor was calculated by matching the inten- 
sities in the two frames of several bright, unsaturated 
field stars, which are assumed not to have Ha emission. 
The accuracy of the scale factor, indicated by the scatter 
shown among the values for different stars, is < 5%. The 
reliability of the method is corroborated by the lack of 
systematic residuals in the stellar images in the net Ha 
frame. 

The next step was the calibration in flux of the 
continuum-subtracted (net) Ha frame. This was done 
by using the observatio ns of spectro photometric stars se- 
lected from the lists of |Oke| (19901. We observed stars 



Feige 34 and bd+33 2642, whose flux is calibrated and 
tabulated in steps of < 2A. This high resolution allows 
us to calculate their flux (or magnitude) through each 
one of the narrow-band filters we used, by integrating 
through their transmission curves. 

2.3. Broadband Imaging 

Broad-band observations of Mrk 35 were carried out in 
January 1999, at the NOT, equipped with ALFOSC. We 



collected CCD images through the B, V, R and / filters. 
The seeing varied between l'.'O and 1'.'5. The log of the 
observations is shown in Table [2j The data reduction 
and calibration is presented in C01a,b. 

2.4. Near Infrared Imaging 

NIR imaging of the sample galaxies was obtained in 
April 2000, at the 4.2-m WHT (ORM, La Palma). We 
used the infrared camera INGRID (Isaac Newton Group 
Red Imaging Device), equipped with a 1024 x 1024 
Hawaii near-IR detector array. The pixel scale is 0'.'242 
pixel^ 1 , and the field of view 4^1 x 4.'1. We took images 
in the J, H and K s bands. The seeing was around 1 
arcsec. A complete log of the observations is provided in 
Table [21 The data reduction is described in C03. 

3. RESULTS 

3.1. Spectroscopic Analysis 

Spectra of Mrk 35 were taken in two almost parallel 
positions (PA. 42° and 38°), as shown in Fig. [T which 
pass through the five brightest SF knots of the galaxy. 
The spatial distribution of the intensity of the bright- 
est emission lines (Ha, H/3, [O III]) has been used as a 
base for defining the regions from where one-dimensional 
spectra of each region were extracted. The different spa- 
tial regions we defined are shown in Fig. [lj the spatial 
length of the slit in Position 1 (SI) was divided into three 
subregions, while Position 2 (S2) was split into two sub- 
regions. 

Figs. [2] and [3] display the integrated spectra of the dif- 
ferent spatial regions, together with the spectrum inte- 
grated over the whole slit length for SI and S2. We 
summarize below the main characteristics of the spectra. 

• Apl includes knot Bl, a markedly red SF region 
(B - V = 0.64; V - I = 0.53) located very close to 
the optical center of the galaxy (see the Ha map in 
Fig. Hp). Absorption wings underneath the Balmer 
emission lines are clearly visible in H/3, H7 and H<5; 
other absorption features such as Ca II A3933, the 
Mg 16 triplet AA5167, 5173, 5184 and Fe I A5335 
are also visible, which indicates the presence of an 
evolved population of stars (> 10 Myr). 

• Ap2 contains knot b2, the optical center of the 
galaxy, which also is a rather red knot (B — V = 
0.63; V — I = 0.70; see Fig. [T|; interestingly, this 
source is a strong emitter in the broad-band frames 

in the NIR it is, indeed, the emission peak — 
but is not an Ha peak. This implies that b2 is 
substantially more evolved. The Ca II A3933, the 
Mg 16 triplet and Fe I A5335 absorption features 
are visible. 

• Ap3 is the sum of the Apl and Ap2. 

• Ap4 corresponds to knot C, a SF region located in 
the galaxy "tail-like" feature. It shows a flat spec- 
trum, with prominent emission lines and no absorp- 
tion features, which is characteristic of a dominant 
OB population. 

• Ap5 corresponds to knot A, the brightest SF region 
of the galaxy. It shows a flat spectrum character- 
ized by strong nebular emission lines, and no vis- 
ible absorption features. The Wolf-Rayet bump is 
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TABLE 2 
Log of the observations 



Date 


Telescope 


Instrument 


Filter/grism 


Exposure time (s) 


Seeing 


Jan. 99 


NOT 2.56m 


ALFOSC 


B 


1800 


l'/55 


Jan. 99 


NOT 2.56m 


ALFOSC 


V 


1200 


to 


Jan. 99 


NOT 2.56m 


ALFOSC 


R 


960 


1'.'44 


Jan. 99 


NOT 2.56m 


ALFOSC 


I 


1200 


l'.'OO 


Mar. 02 


NOT 2.56m 


ALFOSC 


6562 (46) 


3600 


l'/20 


Mar. 02 


NOT 2.56m 


ALFOSC 


6353 (45) 


3600 


1('18 


Apr. 00 


WHT 4.2m 


INGRID 


J 


800 


l'.'OO 


Apr. 00 


WHT 4.2m 


INGRID 


H 


720 


1'.'02 


Apr. 00 


WHT 4.2m 


INGRID 


K s 


1440 


0'.'90 


Feb. 02 


WHT 4.2m 


ISIS 


R300B 


5400 (SI) 


l'/20 


Feb. 02 


WHT 4.2m 


ISIS 


R300B 


3600 (S2) 


1'.'20 



Note. — NOT = Nordic Optical Telescope; ALFOSC = Andalucia Faint Object 
Spectrograph and Camera; WHT = William Herschel Telescope; INGRID = Isaac 
Newton Group Red Imaging Device; ISIS = Intermediate dispersion Spectrograph 
and Imaging System. 



cle arly di scern ible, as already repo rted by |Steel et 



al.|( |1996[ ) and |Huang et al.| ( |1999 l. 



• Ap6 corresponds to knot D, also detached from the 
main body of the galaxy. It too has a flat spectrum, 
without absorption features. 

3.1.1. Reddening Corrected Line Fluxes 

Fluxes and equivalent widths of the emission lines were 
measured using the Gaussian profile fitting option in the 
iraf task SPLOT (a direct integration of the flux for each 
line gave virtually identical results). 

It is known that these measurements are an underesti- 
mate of the real flux of the Balmcr lines, because of the 
underlying absorption component. To corr ect for the un - 



derlying stellar absorption, some authors (McCall, Ryb- 



ski & Shields 1985 Skillman & Kennicutt 1993) adopt 



a constant equivalent width (1.5 — 2 A) for all of the 
hydrogen absorption lines. However, the actual value of 



pends on the age of the star formation burst flDfaz 1988 


Cananzi, Augarde, & Lequeux|1993 


Olofsson 1995 


, anc 



stars. 

The correction from underlying stellar absorption was 
done in two different ways, as explained below. 

For the Balmer lines showing clear absorption wings, 
we fitted simultaneously an absorption and an emission 
component, using the deblending option in SPLOT. 

In the other cases, we proceeded as follows. We first 
adopt an initial estimate for the absorption equivalent 
width, Wabs, correct the measured fluxes, and determine 
the extinction coefficient at A = 4861 A, C(H/3), through 
a least-square fit to the Balmer decrement given by the 
equation: 



F(A) 
F(Hf3) 



J(H/3) 



x 1() C(H/3)x/(A) 



(1) 



where p^p^ is the line flux corrected for absorption and 



normalized to H/3; j^ m is the theor etical value for case 

B recombination, from |Brocklehurst (19711, and /(A) is 
the reddening curve normalized to H/3 which we took 
from|Whitford|dl958l. 



We then vary the value of M^abs, until we find the one 
that provides the best match (e.g. the minimum scat- 
ter in the above relation) between the corrected and the 
theoretical line ratios. This is done separately for each 
knot. 

Reddening-corrected intensity ratios and equivalent 
widths for the different spatial regions, in the two ob- 
served positions, are listed in Tables [3] and [4] 

3.1.2. Physical Conditions and Chemical Abundances 

Table [5] lists the most relevant line ratios, the physi- 
cal parameters and oxygen abundances derived for each 
aperture. 

The excitation ratio [O III]/H/3 is > 2 in all the knots, 
and it is highest in knot A, the brightest star-forming 
region, where it reaches ~ 3.6. [N II]/Hq drops where 
[O III] /H/3 increases, as expected in regions of recent 
star- formation. 

[O I] A6300 is detected in all the apertures, but the 
fluxes measured are quite small. Photoionization by stars 
is the dominant mechanism. Although the presence of 
shocks cannot be completely ruled out, the small values 
of [S II] /Ha and [N II] /Ha indicate that shocks do not 
play a significant role. As expected, in the classical di- 
agrams [N II]/Ha vs [O III]/H/3, [S II]/Ha vs [O III]H/3 
and [O I] /Ha vs [O III] /H/3, all the star- forming knots 
are located within the locus of H II regions. The electron 
density of the ionized gas was determined from the [S II] 
A6717/A6731 ratio. W e used the task temd en, based 
on the fivel program (Shaw & Dufour 19951, which is 
included in the IRAF package NEBULAR. For knot A we 
found a considerably high density (17 cm" 3 ), in good 
agreem ent with the value reported in |Izotov fc Thuan] 
(20041. Knots B, c and D have lower gas densities (< 100 



cm" J ), characteristic of massive H II regions. 

A precise determination of the ionized gas tempera- 
ture (T e ) requires an accurate measurement of the flux 
of [O III] A4363; though we detected this line in five out 
of the six apertures, the uncertainties in the flux mea- 
surements — due to the weakness of the line as well as 
to the strong absorption wings of the nearby H7 line — 
make the derived values unreliable. Only in Ap5 and S2 
the uncertainties in the flux measurements are less than 
25%, and we could compute T c from the line flux ratio 
[O III] A4363/(A4959 + A5007). 
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TABLE 3 

Reddening corrected line intensity ratios 



Line 


on 


Apl 




Ap2 




Ap3 










( A ) 






-W x 








r\ 




-Wx 




3727 


fO III 


2.058 ± 0.043 


92.86 ± 12.74 


2.404 ± 0.050 


78.47 ± 8.53 


2.230 ± 0.024 


72. 


.00 


2 


.20 


3798 


H10 


0.046 ± 0.037 


1 1 A I 

1.14 ± 


u.yo 
















3835 


H9 


0.059 ±0.019 


l.oU it 


U. |Z 


0.044 ± 0.012 


1.08 ± 0.35 


0.048 ± 0.010 


1. 


25 





.30 


3869 


[Ne III] 


0.243 ±0.012 


O.OO it 


U.oo 


0.216 ± 0.012 


3.93 ± 0.27 


0.227 ± 0.005 


4. 


.32 





.12 


3889 


He I + H8 


0.169 ±0.010 


^ c;Q -t 
O.OO it 


U.oz 


0.152 ± 0.011 


4.19 ± 0.42 


0.160 ± 0.007 


4, 


62 





.28 


3968 


Ne III + H7 


0.141 ±0.008 


A 01 -L- 
4.Z1 it 


U.v51 


0.128 ± 0.008 


3.20 ± 0.25 


0.131 ± 0.005 


3. 


,11 





.17 


4026 


He I 


0.031 ±0.018 


U.D ( ± 


n a i 
U.41 
















4069 


SII 


0.025 ±0.011 


U.OO it 


















4101 


H<5 


0.266 ±0.010 


K (17 -t 

o.y i it 




0.297 ±0.011 


5.94 ±0.23 


0.265 ±0.006 


5. 


17 





.13 


4340 


H 7 


0.465 ± 0.005 


14 10-1- 
LHl.LKJ It 


n 99 
u.zz 


0.484 ± 0.008 


11.98 ±0.21 


0.468 ± 0.005 


12. 


.38 





.16 


4363 


[O III] 


0.039 ± 0.011 


1.15 ± 


0.35 


0.032 ± 0.012 


0.75 ±0.31 


0.028 ±0.007 


0. 


.68 





.18 


4471 


He I 


0.064 ± 0.006 


1.77 ± 


0.20 
















4861 


H/3 


1.000 ±0.005 


34.32 ± 


0.15 


1.000 ±0.005 


29.56 ±0.17 


1.000 ± 0.004 


31. 


.47 





.12 


4959 


[O III] 




35.08 ± 


0.22 


0.894 ± 0.005 


25.66 ±0.15 


0.936 ±0.004 


28. 


.16 



u 


.14 


5007 


[O III] 


3.021 ±0.012 


103.10 ± 


U.D1 


2.540 ±0.011 


74.90 ± 0.31 


2.708 ±0.008 


81. 


.65 





.24 


5200 


[NI] 


0.020 ±0.004 


0.80 ± 


0.17 


0.023 ±0.003 


0.78 ±0.12 


0.022 ±0.003 


0. 


.78 





.10 


OO 1 O 


H» T 


0.127 ±0.003 


6.19 ± 


0.22 


0.118 ±0.003 


4.63 ±0.14 


0.123 ±0.002 


5. 


.23 





1 1 


6300 


O I 


0.046 ± 0.004 


2.65 ± 


0.26 


0.061 ±0.004 


2.85 ±0.20 


0.054 ±0.002 


2. 


.66 





.11 


6312 


S III 


0.016 ±0.003 


0.91 ± 


0.20 






0.014 ±0.002 


0. 


68 





.13 


6363 


OI 


0.013 ±0.002 


0.75 ± 


0.14 


0.021 ± 0.004 


1.00 ±0.20 


0.016 ±0.002 


0. 


.78 





.09 


6548 


[N II] 


0.092 ± 0.003 


5.37 ± 


0.28 


0.110 ±0.003 


5.26 ±0.18 


0.101 ± 0.002 


5. 


21 





.13 


6563 


Ha 


2.856 ±0.017 


165.40 ± 


1.96 


2.874 ±0.038 140.40 ± 1.16 


2.862 ± 0.012 146. 


.80 


n 


.85 


6584 


[N II] 


0.292 ± 0.003 


17.30 ± 


0.34 


0.325 ±0.005 


15.78 ±0.26 


0.308 ± 0.003 


16. 


.19 





.20 


6678 


He I 


0.032 ± 0.005 


2.05 ± 


0.36 


0.037 ±0.005 


1.88 ±0.27 


0.032 ± 0.003 


1. 


.75 





.14 


/ 1 I 


re TTi 
fc> 11] 


0.264 ± 0.004 


16.91 ± 


0.37 


0.341 ± 0.006 


17.20 ±0.32 


0.303 ±0.002 


16. 


,18 





. lo 


6731 


[S II] 


0.208 ±0.004 


12.94 ± 


0.34 


0.261 ±0.005 


13.33 ±0.30 


0.234 ±0.002 


12. 


73 




.14 


C(H/3) 




0.050 ± 0.005 




0.064 ±0.016 


0.067 ±0 


.004 






W(HaO„K,(A) 






6.0 




0.0 




5.5 










wens) u (A\ 






3.8 




4.2 




3.5 










VV ^J^JabsV^J 






2.9 




3.5 




3.0 










W(±idJ abs (AJ 






4.2 




3.5 




4.5 














25.2 ±0.4 




27.1 ± 1.1 


55.6 ±0.6 








Line 


Ion 


Ap4 




SI 














(A) 




Fx 


-W X 




Fx 


-Wx 












3727 


fO III 

L w J 






















3798 


H10 






















3835 


H9 








0.039 ± 0.007 


1.04 ± 0.21 












3869 


[Ne III] 


0.212 ± 0.020 


4.81 ± 


U.D4 


0.203 ± 0.005 


4.06 ±0.12 












3889 


He I + H8 


0.188 ±0.021 


6.13 ± 


u.y / 


0.145 ± 0.005 


4.31 ± 0.18 












3968 


Ne III + H7 


0.166 ± 0.019 


4.68 ± 


U.DD 


0.125 ±0.004 


3.48 ±0.12 












4026 


He I 






















4069 


S II 






















4101 


US 


0.343 ± 0.024 


6.81 ± 


U.OU 


0.262 ± 0.004 


5.30 ± 0.09 












4340 


H 7 


0.524 ±0.018 


15.65 ± 


u. / y 


0.471 ± 0.003 


13.41 ± 0.12 












4363 


[O III] 








0.019 ±0.005 


0.47 ± 0.13 












4471 


He I 


0.037 ± 0.012 


0.95 ± 


0.35 


0.049 ± 0.004 


1.19 ±0.11 












4861 


H/3 


1.000 ±0.011 


40.34 ± 


0.74 


1.000 ±0.003 


34.62 ± 0.24 












4959 


[O III] 


0.823 ± 0.009 


29.79 ± 


0.37 


0.878 ± 0.003 


26.88 ±0.10 












5007 


[O III] 


2.399 ± 0.022 


85.79 ± 


U.o4 


2.552 ±0.007 


79.29 ±0.26 












5200 


[NI] 


0.034 ±0.010 


1.44 ± 


0.44 


0.024 ± 0.002 


0.86 ± 0.08 












Oo 1 O 


no ± 


0.119 ±0.007 


6.06 ± 


0.55 


0.121 ± 0.002 


5.02 ±0.07 












6300 


O I 


0.078 ± 0.010 


4.44 ± 


0.73 


0.058 ±0.001 


2.77 ±0.08 












6312 


S III 








0.015 ±0.002 


0.75 ±0.12 












6363 


OI 


0.031 ± 0.010 


1.76 ± 


0.63 


0.018 ±0.001 


0.87 ±0.07 












6548 


[N II] 


0.111 ± 0.009 


6.18 ± 


0.50 


0.104 ±0.002 


5.09 ±0.12 












6563 


Ha 


2.925 ± 0.140 


162.90 ± 


3.04 


2.866 ±0.011 141.30 ±0.73 












6584 


[N II] 


0.327 ±0.017 


17.99 ± 


0.57 


0.314 ±0.002 


15.61 ± 0.16 












6678 


He I 


0.028 ± 0.007 


1.66 ± 


0.50 


0.031 ± 0.002 


1.62 ±0.09 












6717 


[S II] 


0.376 ±0.020 


22.97 ± 


0.78 


0.325 ±0.002 


16.97 ±0.20 












6731 


[S II] 


0.277 ± 0.015 


18.04 ± 


0.71 


0.245 ± 0.002 


12.76 ±0.17 












C(H/3) 




0.035 ± 0.060 




0.001 ± 0.003 












W(Ha) abs (A) 






0.0 




2.5 














W(H/3) abs (A) 






0.0 




2.5 














W(H 7 ) abs (A) 






0.0 




2.3 














W(H5) abs (A) 






5.0 




4.5 














F(H/9) 




8.1 ± 1.2 




57.8 ±0.5 













Note. — Reddening-corrected line intensities, normalized to H/3 = 1, for the measured apertures from the long slit spectrum 
in position 1. Balmer lines are corrected from underlying stellar absorption. The reddening coefficient, C(H/3), the value of the 
absorption correction in the Balmer lines, W abs , and the corrected H/3 flux, .F(H/3) (xl0~ 15 ergs cm -2 s" 1 ). are also listed. 
[O II] data for Ap4 and SI are affected by "ghost emission" from knot Bl and are omitted (see text in § |3.1.2|for details). 
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Fig. 1. — a) top-left: Contour plot of the continuum subtracted Ha image of Mrk 35. The two slit positions are plotted, and the 
different subrcgions we selected from the spectra are marked in bold line, b) top-right: Contour plot of the continuum subtracted Ha image 
of Mrk 35 with labels identifying the individual star-forming knots, c) bottom-left: B-band contours overlaid on the Ha gray-scale map. d) 
bottom-right: i\-band map. Knot b2 is the optical center of the galaxy, but is not a Ha peak. Image orientation is north up and east left. 
Axis units arc arcscconds. 
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Fig. 2. — Spectra of the different subregions defined in slit position Sf, and its total integrated spectrum SI. The label indicates which 
knot each subregion includes. The inset is an enlargement of the spectrum around the H5 line, to enhance the visibility of its absorption 
wings. 
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Fig. 3. — Spectra of the different regions defined in slit position S2, and the total integrated spectrum S2. The inset shows the blue 
Wolf-Rayet bump. 
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TABLE 4 

Reddening corrected line intensity ratios (Position 2) 



Line 


Ion 


Ap5 




Ap6 








S2 




(A) 




r A 




Fx 






''A 




Fx 






3727 


[O II] 


2.023 ±0.032 


117.10 ±4.19 


















3798 


H12 






0.165 ±0 


119 


5 


11 ± 


3.59 








3798 


Hll 


0.085 ±0.013 


5.52 ± 1.43 












0.059 ± 


012 


3.05 ±0.80 


3798 


H10 


0.071 ± 0.008 


3.33 ± 0.48 












0.071 ± 


007 


3.24 ±0.47 


3835 


H9 


0.081 ± 0.004 


4.72 ±0.36 


0.121 ± 


052 


5 


33 ± 


3.27 


0.080 ± 


004 


4.53 ±0.31 


3869 


[Ne III] 


0.274 ± 0.004 


13.77 ±0.22 


0.240 ±0 


028 


8 


21 ± 


1.45 


0.267 ± 


005 


13.23 ±0.17 


3889 


Hel + H8 


0.190 ±0.003 


11.88 ±0.33 


0.189 ±0 


028 


8 


35 ± 


1.77 


0.187 ± 


004 


11.59 ±0.28 


3968 


Nelll + H7 


0.229 ± 0.004 


14.10 ±0.28 


0.197 ± 


022 


7 


73 ± 


1.05 


0.220 ±0 


004 


13.02 ±0.28 


4026 


Hcl 


0.022 ±0.003 


1.17 ±0.19 


0.053 ±0 


047 


1 


57 ± 


1.58 


0.023 ± 


003 


1.14 ±0.16 


4069 


SII 


0.019 ±0.004 


1.01 ±0.22 












0.018 ±0 


003 


0.89 ±0.17 


4101 


H<5 


0.297 ± 0.004 


17.83 ± 0.34 


0.365 ± 


041 


13 


94 ± 


3.18 


0.299 ± 


005 


17.28 ±0.35 


4340 


H 7 


0.489 ± 0.004 


36.25 ±0.50 


0.594 ± 


028 


23 


10 ± 


1.29 


0.494 ± 


005 


34.67 ± 0.43 


4363 


[O III] 


0.025 ±0.002 


1.65 ±0.12 


0.064 ± 


031 


2 


25 ± 


1.30 


0.030 ± 


002 


1.93 ±0.15 


4471 


He I 


0.045 ± 0.001 


3.05 ±0.12 


0.046 ± 


012 


1 


82 ± 


0.58 


0.045 ± 


002 


2.95 ±0.11 


4861 


H/3 


1.000 ±0.002 


90.83 ±0.40 


1.000 ±0 


010 


54 


18 ± 


0.84 


1.000 ±0 


002 


87.46 ± 0.32 


4959 


[O III] 


1.222 ±0.003 


102.50 ±0.36 


0.890 ±0 


011 


42 


73 ± 


0.83 


1.192 ±0 


003 


96.39 ±0.31 


5007 


[O III] 


3.597 ± 0.009 


272.40 ± 0.92 


2.603 ±0 


032 


120 


10 ± 


1.68 


3.506 ±0 


010 


247.70 ± 0.82 


conn 


TNT Tl 

[IN 1] 


0.013 ± 0.001 


1.35 ±0.12 


n nQn -i_ n 


008 


1 


/ Z ± 


U.oo 


U.Ulo ± U 


UU1 


l.oo ± U.lo 


5411 


He II 






















5876 


He I 


0.128 ±0.002 


16.44 ± 0.26 


0.119 ±0 


010 


7 


80 ± 


0.56 


0.126 ±0 


002 


14.84 ± 0.22 


6300 


OI 


0.037 ±0.001 


5.78 ± 0.24 


0.068 ±0 


010 


5 


25 ± 


0.78 


0.039 ± 


001 


5.62 ±0.21 


6312 


S III 


0.014 ± 0.001 


2.06 ±0.14 












0.013 ±0 


001 


1.85 ±0.13 


6363 


O I 


0.012 ±0.001 


1.82 ±0.15 


0.023 ± 


007 


1 


74 ± 


0.55 


0.013 ±0 


001 


1.78 ±0.12 


6548 


[N II] 


0.071 ± 0.002 


9.28 ±0.22 


0.102 ±0 


012 


6 


84 ± 


0.62 


0.078 ± 


002 


10.26 ±0.23 


6563 


Ha 


2.885 ±0.053 


344.00 ± 2.58 


2.998 ±0 


301 


220 


80 ± 


8.17 


2.900 ± 


072 


306.40 ± 1.57 


6584 


[N II] 


0.221 ±0.004 


31.15 ±0.39 


0.317 ±0 


033 


22 


84 ± 


0.98 


0.229 ±0 


006 


30.01 ±0.34 


6678 


He I 


0.034 ±0.001 


5.66 ±0.19 


0.041 ± 


010 


3 


51 ± 


1.07 


0.035 ± 


001 


5.27 ± 0.17 


6717 


[S II] 


0.201 ± 0.004 


32.07 ±0.45 


0.348 ± 


038 


29 


65 ± 


1.63 


0.211 ± 


006 


31.43 ±0.41 


6731 


[SII] 


0.160 ±0.003 


25.88 ±0.33 


0.261 ± 


029 


22 


81 ± 


1.51 


0.168 ±0 


005 


25.65 ±0.32 


C(H/3) 




0.079±0.023 


0.035±0.130 




0.058±0.032 


W(Ha) abs (A) 













6.0 













W(H/3) abs (l) 













0.0 













W(H 7 ) abs (A) 













0.0 













W(H5) abs (A) 













0.0 













F(H/3) 




110.4±6.7 




6.8±2.3 






112.8 ±9.3 



Note. — Reddening-corrected line intensities, normalized to H/3 = 1, for the measured apertures from the long slit spectrum, 
position 2. Balmer lines are corrected from underlying stellar absorption. The reddening coefficient, C(H/3), the value of the 
absorption correction in the different Balmer lines, W abs , and the corrected H/3 flux, F(H/3) (xl0 — 15 ergs crrT_^_s_~ 1 ) , are also 
listed. [O II] data for Ap6 and S2 are affected by "ghost emission" from knot A and are omitted (see text in §|3.1.2|for details). 



An alternative method to esti mate T e is prov ided by 
the empirical relation proposed by Pilyugin (2001 1, which 
relies on the flux of [O II] A3727. 

The [O II] A3727 fluxes we measured in knots C and D 
are abnormally high when compared to their H/3 fluxes 
and give unrealistic values for T c . This problem was 
found to be caused by strong contamination of the [O II] 
line in knots C and D by "ghost [O II] emission" , originat- 
ing from knots Bl and A. This ghost, whose amplitude is 
about one fourth that of the parent [O II] line, is visible 
in all the spectra of the other objects observed in the 
same run. Oddly enough, no ghost is produced by any 
other emission lines. 

Abundances have been estimated from [O II] A3727 
(Pilyugin 2001|) and from the empir ical N2 calibrator, 



log(|JTTOBS84/Ha) introduced by 



ct al. 



et al. 



( 1994 1 ; we used the linear fits der ived in Denicolo 



Storchi- Bcrgmann 



(2002) and Pettini & Pagel (20041. 



3.1.3. Gaseous Kinematics 

The kinematics of the gas was determined by fitting 
a Gaussian to the strongest emission lines (Ha, H/3, 
[O III]A4959 and [O III]A5007). In the regions between 



or outside the knots where the gas emission is fainter, 
we averaged the spectrum by 3 to 5 pixels in the spa- 
tial direction to increase the signal-to-noise ratio. At 
each spatial position, we averaged together the veloc- 
ity values of the measured emission lines (after rejecting 
outliers) , and took their scatter as an estimate of the as- 
sociated uncertainty. A zero point offset on the velocity 
axis was applied to correct for small shifts in the wave- 
length scale, derived by measuring the position of a few 
bright sky-lines. 

The velocity profiles in the two slit positions are very 
similar (Figures [i] and [51. There seems to be an over- 
all, regular rotation in the outer parts of the starburst 
region, with amplitude of about 100 km s , on which 
a seemingly counter-rotating component, with a size of 
about 10", is superposed. This decoupled region is lo- 
cated southwest of the continuum and Ha peaks, and 
can be identified with the depression in gas emission be- 
tween the central regions and the two southwest knots, 
C and D (see Fig. [l| . 

3.1.4. Emission Line Contribution to the Broad-Band filters 
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TABLE 5 

Physical parameters and oxygen abundances for the apertures measured in Mrk 35. 





Apl 


Ap2 




Ap4 


o ± 


An^ 

/\po 








(knot Bl) 


(knot b2) 


(knots B1+B2) 


(knot c) 




(knot a) 


(knot d) 




logQO III]/H/3) 


0.48 


0.40 


0.43 


0.38 


0.41 


0.56 


0.42 


0.54 


logQN II] /Ha) 


-0.99 


-0.95 


-0.96 


-0.95 


-0.96 


-1.12 


-0.98 


-1.10 


log([S II]/Ha) 


-0.78 


-0.68 


-0.73 


-0.65 


-0.70 


-0.90 


-0.69 


-0.88 


log([0 I]/Ha) 


-1.79 


-1.67 


-1.72 


-1.57 


-1.69 


-1.89 


-1.64 


-1.87 


N e (cm" 3 ) 


150 


110 


120 


< 100 


100 


170 


< 100 


170 


Te (K) 












10250 




10770 


T e (K) 


8700 


8800 


8700 






9200 






12 + log(0/H) 


8.46 


8.43 


8.45 






8.43 






12 + log(0/H) 


8.40 


8.43 


8.41 


8.43 


8.42 


8.31 


8.42 


8.32 


12 + log(0/H) 


8.43 


8.36 


8.35 


8.36 


8.35 


8.26 


8.36 


8.27 



Not e. — Line 6: T e derived from th e [O III] A 4363/(A4959 
(|2001k; line 8: abundance derive d from lPilyugin (2001): line 9: 
"TO: c 



line 10: abundance derived fromlPettim & Pagel]T[2004l. 



-I- A5007) ratio; line 7: T fi derived from Pilyugin 
abundance derived from |Denicol6 et al.| ( |2UU2} ; 
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Fig. 4. — Velocity profile for Position 1. The Ha flux profile (in 10 17 ergs cm 2 s 1 arcsec 2 units) is also plotted to help relate 
kincmatical and morphological features. Coordinate R = was set at the center of the peak of the continuum emission, b2. 



In objects with significant gas emission, such as BCDs, 
emission lines can contribute significantly to the flux 
through a broad-band filter, and thus affect broad-band 
colors; the exact amount of such contribution in a given 
band depends on the intensity of the emission lines and 
on their location under the filter transmission curve. 
Most of the evolutionary synthesis models (e.g. STAR- 
BURST99 ) do not take into account the emission line 
contribution to the broad-band filters, and therefore it 



is essential that this contribution be measured and sub- 
tracted out before computing the broad-band colors of 
the knots. 

To do this we proceeded as follows. At each spatial 
position, and for each filter, we computed the integral 
of the curve obtained by multiplying the number of pho- 
tons emitted per unit wavelength interval in the observed 
spectrum by the filter transmission curve. The same cal- 
culation was repeated on the continuum, modeled by fit- 
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Fig. 5. — Velocity profile for Position 2. Coordinate R = was set at the center of knot A. 



ting a polynomial to the spectrum (we found that or- 
ders of about 20 gave the best results) , using an iterative 
sigma-clipping algorithm that rejects all datapoints more 
than about 1 sigma off the fit. 

This was straightforward for the B and V filters, as 
their cutoff values fall inside the wavelength range of our 
spectra. As for R, whose transmission curve extends to 
longer wavelengths than our spectra (which terminate at 
A ~ 6920), we computed R' magnitudes, which we define 
as the magnitude for a R' filter whose transmission curve 
is equal to that of the standard R filter up to A ~ 6920, 
and zero afterward. 

By computing the ratio between the two integral val- 
ues and transforming into magnitudes, we could compute 
the magnitude differences B ^ s — -Bcont > K>bs — Kont an d 
R' ohs — i?cont; as weu as the offsets in color (B — V) b s — 
{B - V) cont and (V - R') ohs - (V - i?')co„t. Note that 
such values are independent of any calibration constant. 

Figs. [6] and [7] display the magnitude difference as well 
as the color offsets, along the two slit positions. From 
these figures it is evident the importance of performing 
a spatial correction: by ignoring the correction or by 
adopting some averaged values, we can obtain misleading 
results. For instance, along position 1 (see Fig. |6|, the 
emission line contribution is < 0.2 mag in V, except for 
knot C (at r ~ 18 arcsec), where it increases to ~ 0.35. 



Along position 2, the emission line contribution is larger 
(with a peak of ~ 0.70 in V in knot a) and shows greater 
variations. 

Unfortunately as we do not have NIR spectra, we could 
not extend this analysis to the NIR domain. 

3.2. Photometric analysis 

3.2.1. Morphological view of Mrk 35 

The B-band contour map and the Ha map of Mrk 35 
are displayed in Fig. [I] The starburst knots are aligned 
along the NE-SW direction, forming a bar-like structure, 
while the brightest knots (a, Bl and b2) form a "heart- 
shaped" structure in the central part of the galaxy. b2 
is the center of the outer isophotes, while knot A is the 
brightest SF region in the galaxy; the peaks of A and B2 
are sep arated by a distance o f 3'.'9 (about 300 pc); in HST 



frames ( Johnson et al.|2004 ) both knots are resolved into 
a number of smaller SF regions, many of them probably 
Super Star Clusters (SSCs). The optical frames suggest 
the presence of a dust lane (crossing the galaxy between 
knot A and knot B, and extending out to the west), which 
is more clearly vis ible in the high resolu tion HST / WFPC 
F606W image by |Malkan et al.| Jl998| ). The dust lane is 
more prominent in the U and is frames, and fade away 
moving toward redder wavelengths; it displays a redder 
color than the LSB component. 
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Fig. 6. — Magnitudes and colors differences between the observed spectrum and the interpolated continuu m for Position 1. The 
discontinuities in the B and B — V plots correspond to the region affected by the " [O II] ghost" , as described in § |3.1.2| 
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Fig. 7. — Magnitudes and colors differences between the observed spectrum and the interpolated continuum for Position 2. 
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A tail-like feature extends in the SW direction, and 
connects with knots C and D, located w 20" (1.5 kpc) 
from the center. The total length of the SF bar-like 
structure is « 70" (5.3 kpc); the equivalent radius of 
the isophote at 27 B-mag arcsec -2 is 56" (4.2 kpc). 

The isophotes twist going from the inner regions, dom- 
inated by the starburst, to the outskirts of the galaxy, 
where they trace the shape of the underlying host, and 
have constant ellipticity for radii > 30". The major axis 
of the outer isophotes has a position angle of 90 degrees, 
and form s an angle of 42 degrees wit h the central chain 



of knots. Sanchez-Portal et al. (2000) reported a P.A. of 



45 degrees for the outer isophotes, while |Steel et al 



(1996) set the major axis of the outer isophotes at P.A. 
ss 58 degrees. In both cases the photometry was not 
deep enough to reach the low surface brightness (LSB) 
galaxy host, and what is claimed to be the "outer parts" 
is actually a region where the starburst emission is still 
significant. This discrepancy stresses once more the need 
for deep, high-quality data to correctly characterize the 
host galaxy in BCDs. 

The morphology of the galaxy in the broad-band opti- 
cal frames basically coincides with that in the infrared, 
though in the optical the inner isophotes are more dis- 
torted, likely due to dust. The spatial position of the five 
stronger sources coincide in all the broad-band frames. 
Knots Bl and b2 are the intensity peaks in the NIR, but 
they are not strong emitter in the blue band; and while 
A is the intensity peak in the bluer bands (U, B, V), it 
is a weak source in the NIR. This fact indicates a highly 
inhomogeneous dust distribution and/or age differences 
among the knots. 

Fig. [TJ; is the gray-scale Ha map of Mrk 35; S-band 
contours are over-plotted. The Ha emission is located in 
the central regions of the galaxy, elongated in the NE- 
SW direction, but slightly off-center to the NW. The un- 
derlying stellar host extends out much farther than the 
starburst region. The peaks of knots A, C and D in the 
continuum (broad-band frames) are also peaks in Ha, 
whereas knot b2, which is a strong emitter in the con- 
tinuum, do not correspond to any of the Ha local max- 
ima. The Ha map reveals intriguing features. Besides 
the eight relatively large SF regions labeled in Fig. [l] 
we found a number of smaller condensations around the 
main body of the galaxy: three knots, together with knot 
E form a " chain- like" structure extending in a north-east 
direction; two other knots are located about 10" north of 
knots C and D; two small condensations are found south 
of knot D 3 . All these knots are continuum emitters, so 
they are stellar clusters. Faint, diffuse extensions emerge 
in the form of plumes and arcs from the main body of 
the galaxy; in particular, a remarkable arc-like structure 
with a diameter about 1 kpc, departs (apparently from 
knot a) towards the West. 

3.2.2. Color maps 

High resolution and deep multi- wavelength broad-band 
observations, from the blue to the NIR, complemented 
by narrow band imaging, are powerful tools to study the 
stellar content of such complex systems as BCDs: i) the 
narrow-band frames allow us to identify the regions of 

3 These knots are faint and not visible in the Ha maps presented 
in Fig. [1] 



the galaxies where active star-formation is taking place 
and, at the same time, the areas free of gas emission; ii) 
in the color maps we can discriminate the different stel- 
lar populations and see what regions are affected by gas 
contamination (for instance in B — V or B — R maps) or 
by dust lanes and patches (optical-NIR color maps). 
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Fig. 8. — B — V color map of Mrk 35; isocontours for the con- 
tinuum subtracted Ha frame are overlaid. North is up and east to 
the left; axis units are arcseconds. 



Fig. [8] displays the B - V map of Mrk 35; we clearly 
distinguish the starburst, elongated in the north-east 
south-west direction and placed atop the underlying ex- 
tended and red elliptical host. The color distribution of 
the starburst region is neither blue nor homogeneous, as 
red patches appears along the chain. This is due to the 
emission line contamination: the strong [O III] AA4959, 
5007 lines account for a large fraction of the flux through 
the V filter, whereas the contribution of the emission lines 
in the B filter is significantly smaller, producing a sub- 
stantial reddening of the B — V color (see Figs, p] and [7]). 
We have over-plotted the Ha emission on the color map; 
notice that the reddest patches inside the starburst coin- 
cide with the stronger gas emitters (namely, A, C and d). 
This result highlights, once again, the need of taking into 
account the contribution of the gas emission and inter- 
nal extinction, when interpreting the broad-band colors 
in starbursts. 

Figs. [9] and [10[ display the B - R and V - K color 
maps ofMrk 35. Because of the large difference in the in- 
terstellar extinction coefficient between blue and infrared 
bands, these maps are useful to recognize and trace dust 
lanes and patches. Unfortunately, the K frame is not as 
deep as the optical, and so we can build reliable color 
maps only for the central regions of the galaxy. A red 
band, crossing the galaxy in the north-south direction, 
between knots A and B, is visible: this is possibly a dust 
lane. The inner regions in these maps are unlike the 
broad-band or the Ha knot distribution, as they display 
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P.A. 
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(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 




B 


20.80 


0.98 


15.04 


22.33 


14.11 


0.72 


74.0 


-0.02 




V 


19.45 


1.01 


14.72 


21.58 


13.41 


0.71 


77.3 


0.03 




R 


19.45 


0.97 


14.94 


21.28 


13.11 


0.70 


77.7 


-0.01 




I 


14.50 


0.83 


15.33 


20.80 


12.63 


0.71 


74.9 


0.01 




J 


12.80 


1.12 


16.70 


20.52 


11.65 


0.70 


58.1 


0.18 




H 


14.02 


1.01 


19.99 


20.00 


11.19 


0.67 


63.2 


0.15 




K 


12.80 


0.88 


20.50 


19.70 


10.89 


0.70 


60.6 


0.03 
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Fig. 9. — B — R color map of Mrk 35. North is up and east to 
the left; axis units are arcseconds. 
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Fig. 10. — V — K color map of Mrk 35. North is up and east to 
the left; axis units are arcseconds. 



a patchy pattern, which probably results from the mix- 
ture of gas and dust. 

3.2.3. The Low Surface Brightness Stellar Component 
Surface photometry of Mrk 35 in the optical and in 



the N IR was presented in COlb, C03 and Caon et al. 
(2005). In these works, the 1-D luminosity profile of 



the L SB component was modeled by a Sersic law ( |Sersic| 
1968 1. Here, we carry out a two-dimensional modeling of 



TABLE 6 

2D Sersic photometric parameters for the LSB host galaxy 



Note. — Column (1): filter; col. (2): equivalent radius of the mask 
covering the central starburst region (arcsec); col. (3): Sersic shape in- 
dex n ; cols (4) and (5): effective radius (arcsec) and surface brightness 
(mag arcsec -2 ); col. (6): total magnitude of the LSB component de- 
rived from the fit (mag); col. (7): axis ratio; col. (8): position angle; 
col. (9): diskiness (negative)/boxiness (positive) parameter. 

the surface brightness distribution of the LSB componen t 
following the method described in Amorin et al. ( 2007 1 . 
The two main advantages of this novel approach are: the 
possibility of masking out the starburst region by follow- 
ing exactly its shape, rather than setting an inner radial 
limit to the 1-D profile fit; and the fact that we can then 
subtract out the best fit model from the galaxy image to 
recover the starburst. 

The fit is done on the outer parts of the galaxy after 
masking out the central starburst and other disturbances 
such as small knots, foreground and background objects, 
etc. We used both Ha images and color maps to measure 
the size of the central starburst region. We start with a 
first mask, and run the program. Discrepant points are 
flagged out and the fit is re-done. The process is iterated 
until the fit converges, within a preset tolerance, to a 
solution that does not depend on the exact values of the 
input parameters. 

The resulting two-dimensional Sersic parameters of 
Mrk 35 in the optical and in the NIR bands are shown in 
Table [6] The model built in this way is finally subtracted 
from the original galaxy frame to recover the young star- 
burst. 



3.2.4. Disentangling the light from the different components 

The flux received from any given point in the galaxy is 
the sum of three components: the low surface brightness 
host, the continuum flux from the young stellar compo- 
nent and the flux in emission lines: 



F u 



LSB 



young 



cmlincs 



(2) 

The photometric analysis allows us to compute Flsb 
by means of the Sersic modeling, as explained in the 
previous section, and thus derive 



-^young -^cmlines F\i 



F 



LSB 



(3) 

spectroscopic analysis allows us to compute 
by modeling the spectral continuum as 



The 

Flsb-I 

described in Sect |3.1.4| and thus derive 



youn. 



ilincs 



LSB 



F v , 



(4) 



Therefore, by combining photometry and spectroscopy 
we are able to derive separately each of these three com- 
ponents in filters B, V, and R (the latter if we agree 
that the difference between the R' and the R band is not 
significant for our purposes.) 
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TABLE 7 
Flux decomposition for 
Mrk 35 SF KNOTS 



TABLE 8 

Ho PHOTOMETRY OF THE INDIVIDUAL KNOTS IN MRK 35 



Component B 


V 


R 




Knot A 






p mlin 


24.4 


44.9 


39.8 


^LSB 


10.5 


10.0 


5.4 


■f young 


65.1 


45.1 


54.8 




Knot Bl 






p mlin 


8.0 


16.0 


15.2 


^LSB 


8.5 


8.2 


4.1 


■fyoung 


83.5 


75.8 


80.7 




Knot C 






-^cmlincs 


15.2 


22.9 


24.8 


^LSB 


39.3 


38.2 


18.4 


■fyoung 


45.5 


38.9 


56.8 




Knot D 






-^cmlincs 


19.2 


29.7 


32.6 


^LSB 


36.9 


34.0 


16.4 


■fyoung 


43.9 


36.3 


51.0 




Knot b2 






-^cmlincs 


7.0 


12.9 


14.0 


^LSB 


3.6 


3.3 


3.9 


■fyoung 


89.4 


83.8 


82.1 


Note. — Values are normalized 
to a total flux F to t = 100. 



We illustrate how we proceed by taking knot A, whose 
diameter is about 6". In A the underlying stellar com- 
ponent, as described by the Sersic model, accounts for 
~ 10% of the total flux in B and in V, and 5% in R. In 
the spectral aperture 6" wide centered on this knot, the 
continuum (that is, Flsb + -Fyoung) accounts for 75% of 
the flux in B, the emission lines 25%. In V we have 55% 
and 45% respectively; in R' 60% and 40%. 

Combining the above data, under the simplifying as- 
sumption that the spectral data are constant within the 
whole A area, it is straightforward to derive that the 
young stellar component, F yonnE , accounts for 65% of 
the total flux in B. 

The same procedure can be applied to the other knots. 
The results are summarized in Table [T] 

3.2.5. Photometry of the starburst 

In order to identify the individual SF knots we used 
the FOCAS package and applied it to the continuum- 
subtracted Ha images. FOCAS looks for local maxima 
and minima in the counts of the pixels in the image; we 
adopted the criteria that, for a knot to be detected, its 
area must be larger than 38 pixels (corresponding to an 
equivalent diameter of ta l'/3 or 100 pc), to ensure that 
the diameter of the knot is larger than the point spread 
function; also the counts of each of its pixels must be 
higher than 2.5 times the standard deviation of the sky 
background. Using these criteria eight knots were iden- 
tified (see Fig.[l]D). 

The total Ha flux of the galaxy is 1.09 x 10 -12 ergs 
cm~ 2 s _1 (after setting a threshold of 3.8 x 10 -18 ergs 
cm~ 2 s _1 , 2.5 times the standard deviation of the sky 
background). With a diameter of 0.64 kpc and a total 
flux of 5.57 x 10~ 13 ergs cm -2 s" 1 , A is the brightest Ha 
knot, emitting half of the total Ha flux of the galaxy. 

Next, we computed the broad-band magnitudes of the 
starburst knots identified in the Ha frame, after con- 
volving the broad-band images to match the PSF of the 



Knot 


Flux(Ha) 


1 T /TT \ 

log L(tia) 


— W(Ha) 


-W(Ha) cor 


A 


556.60 


40.21 


450 


572 


Bl 


44.54 


39.11 


157 


178 


C 


23.98 


38.84 


309 


687 


D 


15.86 


38.66 


361 


714 


E 


0.27 


36.90 


151 




F 


6.60 


38.28 


166 


244 


G 


2.02 


37.77 


72 


95 


H 


6.04 


38.24 


102 


143 


B2 






124 


140 



Note. — Values are corrected from interstellar extinction 
and from [N II] emission. The equivalent widths are shown 
before and after the correction from the contribution of the 
LSB component to the continuum. Ha fluxes are in 10 -15 
ergs cm -2 s -1 units; Ha luminosities in erg sec -1 units. 
W(Ha) is in A. 

image with the worst seeing. Results of the photometry 
are presented in Tables[8]and [9j Table [8] displays the Ha 
flux, luminosity and equivalent width of the knots. The 
central knot B2, not detected by FOCAS, is also included 
(fluxes are computed within a circular aperture 1'.'5 in 
radius). The data have been corrected from interstellar 
extinction and [N II] contribution using our spectroscopic 
information. 

Table |9] lists the broad-band colors of the knots cor- 
rected for Galactic extinction, together with the "true" 
starburst colors, computed by subtracting th e cont ribu- 
tion of the underlying host galaxy (see Sect 3.2.3) and 
of the flux due to em ission lines (for the B, V and R 
filters, see Sect 3.2.4), and by correcting for interstel- 
lar reddening. For those knots whose position matches 
a spectroscopic aperture (namely A, Bl, b2, C and d) 
we correct from reddening and emission lines using the 
values derived for the corresponding aperture; in the re- 
maining cases we use the values derived for the nearest 
subregion in the spectrum. 

3.3. Disentangling the Stellar Populations 
3.3.1. The Starburst 

In order to constrain the properties of the SF knots 
identified in Mrk 35, we ad opted the Starburst 99 evo- 
lutionary synthesis models ( Leitherer et al.|1999 l, which 
produce a comprehensive set of observables of galax- 
ies with active star formation. We use models with 
z = 0.004 and z — 0.008, which bracket the metal- 
licity derived from the emission line fluxes (which is a 
reasonable approximation to the metallicity of a young 
population), and look for the scenario that best repro- 
duces the collection of observables available for every 
knot: W(Ha), the optical and NIR colors, as well as 
W(Hf3) for those knots with spectroscopic information. 

In all cases, we can reproduce the Ha equivalent width 
with an instantaneous burst of star formation and the 
Salpeter IMF with an upper mass limit of 100 Mq. Knot 
A, C and D are the youngest; their observables are con- 
sistent with ages of about 5 Myr. The central knots, Bl 
and b2, arc slightly older, with ages of about 7 Myr. 

In Fig. 11 our results are compared with the mod- 
els. The equivalent width of the Ha line is plotted ver- 
sus B — V . All the knots have Ha equivalent widths in 
good agreement with model predictions, with ages rang- 
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TABLE 9 

Observed broad-band photometry of the individual knots detected 

in Mrk 35. 



Knot Area B B - V V - R V - I V - J J - H H — K 



Total magnitudes and colors 



A 


25.7 


15.43 


0.66 


0.01 


0.15 


0.92 


0.41 


0.49 


Bl 


4.6 


17.21 


0.64 


0.12 


0.53 


1.31 


0.42 


0.46 


C 


8.4 


18.61 


0.61 


0.09 


0.35 


0.98 


0.50 


0.31 


D 


5.4 


19.30 


0.67 


0.07 


0.30 


1.00 


0.49 


0.43 


E 


2.4 


22.33 


0.67 


0.26 


0.78 








F 


4.9 


19.38 


0.49 


0.20 


0.54 


1.27 


0.38 


0.42 


G 


2.7 


19.09 


0.50 


0.23 


0.62 


1.33 


0.47 


0.35 


H 


5.6 


18.61 


0.49 


0.25 


0.63 


0.95 


0.49 


0.39 


B2 


7.1 


16.89 


0.63 


0.14 


0.70 


1.43 


0.49 


0.42 



Magnitudes and colors for the young component 



A 


15.68 


0.22 


0.18 


-0.01 


0.66 


0.45 


0.55 


Bl 


17.27 


0.50 


0.17 


0.46 


1.15 


0.46 


0.49 


C 


19.37 


0.41 


0.48 


0.04 


0.62 


0.47 


0.44 


D 


20.10 


0.43 


0.41 


-0.02 


0.70 


0.35 


0.66 


E 


24.40 


0.75 




0.85 








F 


19.66 


0.23 


0.42 


0.40 


0.78 


0.47 


0.55 


G 


18.25 


0.35 


0.36 


0.54 


1.02 


0.60 


0.43 


H 


18.87 


0.22 


0.47 


0.54 


0.62 


0.62 


0.46 


B2 


16.83 


0.52 


0.08 


0.63 


1.28 


0.52 


0.44 



Note. — Upper rows: broad-band photometry for the star-forming regions 
detected in Mrk 35. Aperture photometry (within a radius = 1'.'5) for knot B2 
is also included. Integrated total magnitudes have been corrected for Galactic 
extinction. Bottom rows: Magnitudes and colors of the young component, 
obtained by subtracting the low surface brightness host in every filter, and the 
contribution of emission lines in B and in B — V and V — R colors. Magni- 
tudes and colors have also been corrected for both Galactic and interstellar 
extinction. 



ing from 4 to 7 Myr, but lie far from the models in the 
B — V axis. 

This is not surprising. Although the data have been 
corre cted for the contribution of emission lines (Sec- 
tion |3.1.4[ ), the contribution of the LSB host (Sec- 
tion |3.2.3p and for interstellar extinction, the various 
assumptions and simplifications we made to be able to 
work out such effects (for instance, that C(H/3) and the 
emission line contribution do not vary spatially within a 
same knot) result in significant uncertainties on the final 
values. Also, in Mrk 35, the presence of a substantial 
amount of dust can also play an important role. 

Fig. [12] displays the nir color-color diagram of the SF 
regions. All the observed knots lie off the region of the 
models. We must keep in mind that NIR colors have not 
been corrected from gas line emission, and that SB 99 
includes only the contribution of the nebular continuum. 
However such a large disagreement is difficult to explain 
only in terms of emission lines. These excesses in the IR 



signature of the presence of hot dust (Vanzi 2002; Vanzi 


& Sauvage 2004 Johnson et al.|2004 


Cresci et al.|2005). 



3.3.2. The Low Surface Brightness Stellar Component 

To constrain the evolutionary status of the host galaxy 
in Mrk 35 we can only use the information derived from 
its colors. Although some absorption features are de- 
tected in the spectra, they cannot be measured with suf- 
ficient accuracy. 

The colors of the underlying galaxy host (de rived from 
the Sersic models) are tabulated in Table 10 They are 
not corrected from internal extinction, because the coef- 
ficients we derived, based on the gas emission lines, may 




-0.2 0.2 0.4 0.6 0.1 

B - V 



Fig. 11. — Logarithm of the equivalent width of the Ha line 
versus (B — V). Solid points are the observables of the the knots 
after applying the corrections from reddening and from the contri- 
bution of emission lines and older stars. We also plot the tracks of 
two instantaneous bursts of z = 0.004 (solid line) and z = 0.008 
(dashed line), both aged from 1 to 30 Myr, and with an IMF with 
a = 2.35 and M up = 100M Q (Leitherer et al. 1999). The numbers 
along the tracks mark the age in Myr. 



not apply to the regions outside the area occupied by the 
starburst. 

We use the spectral energy distribution (SED) derived 
from the Sersic model to estimate the age of the under- 
lying old stellar compon ent. We employ models derived 
from the GALEV code (iBicker et al. 1120041 iWeilbacher 
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Fig. 12. — J — H vs H — K color-color diagram. Solid points 
represent the knots parameters after applying the corrections from 
reddening and older stars. We also plot the tracks of two instanta- 
neous bursts of z = 0.004 (solid line) and z = 0.008 (dashed line), 
both aged from 1 to 30 Myr, and with an IMF with a = 2.35 and 
M up = 100 Mq (Leitherer et al. 1999). The numbers along the 
tracks mark the age in Myr. 



TABLE 10 

Colors of the LSB stellar component derived 

BY FITTING A SeRSIC MODEL 

B — V V - R V - I V - J J - H H - K 



0.70 



0.30 



0.78 



1.76 



0.46 



0.30 




SSP, 1.1 Gyr, Z=0.008 
t.1 Gyr, 2.6 Gyr, Z=0.008 
Sb, 8.7 Gyr, cc 
observed values f 

1.2 1.4 
Wavelength [um] 

Fig. 13. — The B to K spectral energy distribution of the ob- 
served LSB component and three typical models. The y-axis shows 
the color relative to V, the x-axis the wavelength. The observa- 
tions are shown as crosses, where the x-error gives the approximate 
width of the filters and the y-axis the uncertainty of the Sersic fit 
in each filter. Three selected models are shown: 1) an SSP with 
Z = 0.008 at an age of 1.1 Gyr; 2) a model with characteristic star- 
formation timescale r = 1 Gyr with Z = 0.008 at an age of 2.6 Gyr; 
3) a chemically consistent model of an Sb galaxy (t = 6.6 Gyr) at 
an age of 8.7 Gyr. 



as described by Anders et al. (20041. Our model grid 
consists of single metallicity models with Z — 0.004 and 
Z = 0.008 and we use single stellar population (SSP) as 
well as models with exponentially decreasing star forma- 
tion rate with a characteristic timescale 1 < r < 12 Gyr. 
As an alternative we also include chemically consistent 
models similar to those of Bicker et al. ( 2004 ) , which are 
divided into different star formation histories to repre- 
sent all galaxy types observed at the present day and use 
exponential timescales from r = 1 Gyr to constant star 
formation rate. 

To look for the best agreement between the measured 
colors and the models we then derive an error weighted 
X 2 value of the observed magnitudes of every time-step of 
each model. Several models give a good fit; all those with 
half solar metallicity give a better fit than those with Z = 
0.004. It is, however, difficult to obtain a good estimate 
of the dominant age of the stellar populations. Despite 
the wide wavelength range covered, the possible solutions 
are dege nera te regarding age and star formation history 
(see Fig. 13 1. While the best matching single metallicity 
model is the SSP with Z = 0.008 at an age of 1.1 Gyr, 
other models with longer star formation timescales and 
only slightly worse x 2 values match at ages of 2.5 to 
7.9 Gyr. Some of the chemically consistent models also 
give very go od fits, the Sb-galaxy modeled after picker 
et al. ( 2004 1 provides the best overall match at an age of 
KT&yr. 

We cannot help concluding that further observations, 
and especially data on spectral indices in the outer part 
of the LSB component, are necessary to determine the 
true dominant age of the underlying stars in the galaxy. 

4. DISCUSSION 

This work is the second of a series of papers devoted 
to a thorough spectrophotometric analysis of a sample 
of nearby BCDs. As we concluded in our previous work 
(C02, focused on the i,E BCD Mrk 370), this kind of 
analysis, although scarce in the literature because of its 
difficulty, is essential to disentangle the stellar popula- 
tions in systems as complex as BCDs and to put strong 
constraints on their SF histories and evolutionary status. 

Here, by combining a wide set of photometric and spec- 
troscopic data we disentangled the two stellar popula- 
tions in the iE BCD Mrk 35: the starburst, which forms 
a bar-like structure extending in the northeast-southwest 
direction, and the underlying population of substantially 
older stars. The current star-formation not only takes 
place in the central k nots (a and b), as pointed out by 



Johnson et al. (20041 — who worked with frames that 
covered only the innermost regions — , but also in the 
two knots placed at the end of the "tail-like" feature (c 
and d), as their Ha emission reveals. 

The brightest knot, A, accounts for half of the Ha emis- 
sion of the galaxy; the age we derived for it is about 5 
Myr, consistent with the ages found in previous works 



(Johnson et al 
detection of W 



2004 Steel et al. 1996 , and with the 



stars ( |Steel et al.||1996| Huang et al 



et al. 2004); all of them use a Salpeter initial mass func- 
tion and Padova isochrones and other input parameters 



19991. Knots C and D show almost identical ages, which 
rules out the hypothesis that the star formation is prop- 
agating along the "tail-like" feature. Knots Bl and B2 
have slightly larger ages, about 7 Myr, but the age differ- 
ence with the brighter knots (~ 2 Myr) is probably too 
small to be consistent with a scenario of self-propagating 
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star- formation. 

The high resolution frames combined with spatially re- 
solved long-slit spectroscopy allowed us to carry out an 
exhaustive analysis of the starburst population. In order 
to characterize the current starburst episode, a consider- 
able effort has been made to disentangle the young stars 
from the ionized gas emission, from dust and from older 
stars. Using deep broad-band frames we modeled the 
LSB host and subtracted it from the total emission; the 
spectra allowed us to estimate the contribution of emis- 
sion lines, and to derive the extinction coefficient, in the 
different spatial regions along the two slit positions (see 
Table [tJ. 

After applying these corrections, we derived colors, 
physical parameters and oxygen abundances for the SF 
knots, and investigated how they vary from one spatial 
location to another. 

A very important aspect, often ignored in previous 
analyses, is the estimation of the contribution of emission 
lines to the total amount of flux measured using standard 
broad-band filters. In this case, using our spectra we cal- 
culated the amount of flux due to emission lines in the 
the B, V and R bands, along both slit positions. We 
found that this contribution is by no means negligible, 
and even more importantly, it shows spatial variations, 
not only from one knot to another, but even within the 
same knot. For instance, the contribution of emission 
lines in knot A (along the slit position 2), is as large as 
0.7 mag; such a huge amount of flux can make the B — V 
starburst colors similar or even redder than those of the 
host (see the B — V color map in Fig.[8|. Therefore, for 
a proper subtraction of the emission lines contribution a 
two-dimensional mapping of the line emission is required. 

Several recent works have questione d the previous be- 



lief that BCDs have no or little dust ([Hunt et al.||2001 
Cairos et aLp 003; Va nzi fc Sauvage||2004| ICabanac et 



aT 2005 1. Mrk 35 shows evidence of dust in its central 



regions. In the optical frames dust is visible as a lane be- 
tween the central knots, A and B, extending to the wes t 
(as already pointed out by C03 and Johnson et al.|2004 1 . 
Our color maps also suggest the existence of a more ex- 
tended dust distribution: red patches are in fact visible 
at the north of knot A, and the pronounced "red" border, 
crossing the galaxy parallel to the central like-bar star- 
burst, may also be due to dust which blocks light from 
the part of the galaxy behind. 

One intriguing aspect is that, although we found ev- 
idence for the presence of dust, the derived reddening 
values in all the apertures are considerably small. From 
the Balmer decrement, we obtain Ay = 0.17 in knot 
A and Ay — 0.11 — 0.13 in knots Bl and B2 (in good 
agreement with the v alues reported in |Steel et al.||1996| . 
Johnson et al.l ( |2~004| ) however, working in NIK and radio 



wavelengths, found considerably higher extinction value 
for knot A, Ay — 8. These discrepancies in the extinc- 
tion measured at different wavelen gths have been also 



Yun et al.|1994{ He2- 



fou nd in other SF dwarfs (V V 114, 

10, |Kobulnicky et al.||1995) : in a dusty environment, m 
which the dust is mixed with the emitting sources, there 
is a trend of increasing extinction from the optical to the 
IR. Therefore, the lower extinction values derived from 
our optical spectroscopy are most probably the result of 
a inhomogeneous dust distribution. 
Indeed, the actual amount of dust (and its distribu- 



tion) is currently a hot topic in BCD galaxy research. 
By studying the properties of the dust in these nearby 
dwarfs we can learn about the dust formation in primor- 
dial environments. Moreover, it is becoming more and 
more clear that an important fraction of the SF activity 
in these gala xies lies buried in dust and has not visible 
counterparts (Thuan et al. 1999 Kobulnicky fc Johnson 



19991 IHunt et al.||2001[ | Vacca et al.||2002| |Johnson et 
il.||2003l |Vanzi fc Sauvage||2004p. Integral Field spectro 
scopicobi 



scopic observations, with which to build a spatial map of 
the extinction coefficient, together with high-resolution 
MIR observations to map the dust distribution and mea- 
sure its emission, would provide important advances in 
the field. 

The galaxy morphology — two major knots in the in- 
ner regions of the galaxy and the "tail-like" features — 
suggests that an interaction (and/or) merger could have 
played a main rol e in the galaxy shaping and evolution. 



Steel et al. ( 1996 1 rejected the merger hypothesis on the 
basis on the velocity of the central knots (almost identi- 
cal) and the regularity of the outer isophotes, and favored 
the hypothesis of intrinsic self-induced star-formation. 
Nevertheless, only small age gradients have been detected 
across the galaxy, eve n with the deeper and better spa- 
tial resolution data in Johnson et al. (20041 and in this 
work. Although we cannot discard that some propaga- 
tion is operating in Mrk 35, it does not see ms to be the 
dominant mechanism. On the other hand, Johnson et 



al. ( 2004 1 favored a merger event as the one igniting the 
star-formation in Mrk 35: they argued that the regu- 
larity of the outer isophotes does not exclude a small- 
scale interaction, and the fact that A and B had similar 
rcdshifts also does not exclude that they were once sep- 
arated systems. They compared Mrk 3 5 with the iE, 



BCD He 2-10, a very well studied object flKobulnicky et 
al.|1995||Johnson et al |2000| ICabanac et al.|2005p , which 
has significant similarities with Mrk 35, and appears to 
be interacting with a massive clo ud of gas that is f alling 
into the main body of the galaxy ( Kobulnicky et al.|1995 

lT~u,~~.-.~ ~i lonnnh ~. 1^.^ \i~ i-u. 



| Johnson et al.|2000"| ) , and speculate that the nature of the 
star-forming trigger in Mrk 35 could be similar. 

Our kincmatical analysis seems to support the merger 
hypothesis. The apparently counter-rotating feature in 
the inner region of the velocity profiles could be the sig- 
nature of a past merger or acquisition event. Integral 
field spectroscopy, which allows the derivation of the ve- 
locity field of the gas, is required to properly interpret the 
rotation pattern. Also, deep spectroscopic observations 
with large telescopes, which permit high signal-to-noise 
observations of absorption features and therefore to trace 
the kinematics of the stars, would be essential. 

5. SUMMARY AND CONCLUSIONS 

Optical and NIR broad-band images, Ha narrow-band 
frames and long-slit spectroscopic data for the i,E BCD 
galaxy Mrk 35 have been analyzed in order to derive the 
properties of the different components of the galaxy, to 
constrain its evolutionary status and to investigate the 
possible mechanisms triggering the actual SF episode. 
From this study we highlight the following results: 

• Two different stellar components are clearly distin- 
guished in the galaxy: the actual starburst and an 
underlying older population. The current star for- 
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mation activity takes place in several knots, aligned 
along the NE-SW direction, in a bar-like structure; 
the brightest knots (a, Bl and B2) are arranged in 
a "heart-shaped" structure, in the central region 
of the galaxy, while knot C and D are located in 
the "tail-like" feature. An extended LSB compo- 
nent, with regular appearance and red colors, sits 
underneath the star-forming area. 

• The brightest knot, A, is a powerful and young 
SF region, with Wolf-Rayet stars, and accounts for 
50% of the total Ha luminosity of the galaxy. 

• For each spatial region in the slit we derived red- 
dening corrected intensity ratios, extinction values, 
physical parameters of the gas and oxygen abun- 
dances. Although we found evidences for the pres- 
ence of dust, the extinction coefficient derived from 
the Balmcr decrements are similar and small for all 
the regions {Ay ranges from to 0.17). The oxy- 
gen abundances do not show significant variations 
from knot to knot. 

• The observed flux in the SF knots is the sum of, 
basically, three components: the light coming from 
the young stars, the contribution of emission lines 
and the emission from the low surface brightness 
underlying population of stars. We presented a 
methodology to disentangle and assess each of these 
components. 

— Using our spectroscopic information we first 
computed the contribution of emission lines 
to the B, V and R filters, along the two slit 
positions. We found that they can strongly 
affect the broad-band magnitudes and col- 
ors, as much as 0.7 mag in the U-band fil- 
ter. This contribution also varies significantly 
from one knot to the next. Thus, only by a 
two-dimensional mapping of the gas emission 
(which can be done by means of Integral Field 
Spectroscopy) can we derive accurate photo- 
metric values of the young star colors. 

— We modeled the host galaxy with a Sersic 
function, and subtracted it from the total 
emission to recover the colors of the young 
SF knots. 



• We compared the final observables with the predic- 
tions of evolutionary synthesis models, and found 
that we can reproduce the colors of the knots with 
an IB of star-formation and the Salpeter Initial 
Mass Function with an upper mass limit of 100 
M©. The ages of the knots are about 5 Myr, with 
B being slightly older. No significant age gradients 
have been found in the starburst region, which 
discards the propagated star formation scenario. 



We found some evidence for the presence of dust in 
the central regions of the galaxy. 



The kinematical analysis shows an overall rotation, 
with a peak-to-peak amplitude of <~ 100 km s _1 . 
However, the curve is irregular in the inner region, 
where there seems to be a counter-rotating compo- 
nent. This could support the idea that the galaxy 
has undergone or is undergoing an interaction or 
accretion process. 
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